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Abstract: Massive gravity can be described by adding to the Einstein-Hilbcrt action a function V 
of metric components. By using the Hamiltonian canonical analysis, we find the most general form 
of V such that five degrees of freedom propagate non perturbatively. The construction is based on a 
set of differential equations for V , that remarkably can be solved in terms of two arbitrary functions. 
Besides recovering the known "Lorentz invariant" massive gravity theory, we find an entirely new class 
of solutions, with healthy features on the phenomenological side, in particular they are weakly coupled 
in the solar system and have a high ultraviolet cutoff A 2 = (m-M*,;) 1 / 2 , where m is the graviton mass 
scale. 
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1 Introduction 

Whether general relativity (GR) is an isolated theory is an interesting question both from the theo- 
retical and phenomcnological side. In particular, the long standing quest for a healthy modification 
of gravity at large distances recently has been target of a renewed interest. The attempt to find 
non-derivative, namely massive, deformations of GR dates back to the work of Fierz and Pauli [1] 
(FP), where it was realized that, around Minkowski space, a generic Lorentz- invariant (LI) massive 
deformation of gravity propagates 6 degrees of freedom (DoF) already at quadratic level, and that the 
sixth mode is a ghost. While one can tune the mass term to get rid of the sixth mode at quadratic 
level, it reappears at nonlinear level or around non-flat backgrounds [2]. Surprisingly, for a long time, 
a complete canonical study to understand the propagating degrees of freedom non perturbatively, in a 
generic deformed version of GR was missing. The gap was filled in [3] where it was shown that five DoF 
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propagate if the deforming potential solves a Monge- Ampere supplemented with an extra equation 1 . 
As matter of fact, infinitely many solutions exist, and the only one discovered in recent years [5-7] is 
just one of them. The peculiarity of this theory is that it exhibits an accidental Lorentz symmetry 
around a single background, usually taken to be the flat one. Conversely, most of the massive gravity 
theories with five DoFs are intrinsically Lorentz violating. As discussed in section 2, this possible 
residual Lorentz symmetry acts just in the pure gravity sector, and is not the same symmetry invoked 
by the Einstein equivalence principle, which is phcnomcnologically required. In fact, it is actually the 
opposite, this symmetry is rather costly from a phenomenological point of view: the Lorentz-invariant 
theory at quadratic level reduces to FP and it shares with FP the failure to reproduce the correct light 
bending, in sharp contrast with GR. The issue does not disappear in the limit of vanishing graviton 
mass, the so called vDVZ discontinuity [8]. A possible way out was proposed by Vainshtein in [9], 
by arguing that the linear approximation cannot be applied in the solar system and nonlinear effects 
restore the correct GR behaviour at short distances from a source, or in the limit of vanishing graviton 
mass. While this mechanism has been verified in a number of models [10], it is important to point 
out that for it to work, the theory has to rely on strong nonlincarities even at the macroscopic scales 
of the solar system where the values of the gravitational potential is small. Quantum corrections also 
seem to be important at macroscopic scales [11]. This makes the theory essentially hard to handle, 
since one can not use perturbation theory even in a two body problem like the motion of a planet 
around the Sun. 

As we show, this problematic situation is avoided if Lorentz invariance in the gravitational sector 
is not imposed. As this concerns only the gravitational sector and the equivalence principle still 
can be preserved, this scenario is not only phcnomcnologically allowed but may also be considered 
interesting. At the quadratic level, a class of such theories that are free from the vDVZ discontinuity 
are available [12-16]. Remarkably many of the solutions of the Monge- Ampere problem arc precisely 
of this type. 

In this paper we give the construction based on the canonical analysis that leads to the general 
massive gravity theories with five propagating DoFs. It turns out that the deforming potential has to 
satisfy a set of differential equations, whose general solutions can be found in terms of two arbitrary 
functions. The result is also important from a phenomenological point of view; indeed, it provides a set 
of new theories [17] which are weakly coupled in the solar system and with a cutoff A2 = (mMpj) 1 / 2 , 
the largest possible in the absence of a fundamental Higgs mechanism for gravity. 

The outline of the paper is the following: in section 2 wc discuss the various formulations of 
massive gravity, stressing the important role of the residual symmetries of the deforming potential V. 
In section 3 a detailed study of the number of propagating DoFs is given for a generic nonderivative 
modification of GR encoded in V; the result is that five DoF are present if V satisfies the Monge- 
Ampere equation supplemented by an extra differential equation in field space. In section 4 we find 
the general solution to these equations, which leads to the most general potential V which propagates 
five DoFs. In section 5 the potential V is constructed explicitly in a number of cases. In section 7 the 
perturbation theory around Minkowski space is considered. Finally, in section 8 we summarize our 
conclusions. 



Sec [4] for a alternative analysis using Kuchar's Hamiltonian formalism. 
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2 Massive GR: From Lorentz Invariance to Rotational Invariance 



To deform GR into a massive theory we need to add nondcrivativc functions of the metric field 
Because any scalar diff- invariant function of the metric is trivial, this can only be done by breaking 
cliffs and thus reducing the invariance properties of the action. The best one can do is to preserve 
a residual set of symmetries in some specific background(s). Indeed, formally this can be described 
through the addition of fixed, external, tensor quantities, with different kind of resulting theories. 2 

• A traditional way was to add an external second metric g llVl so that combining g and g one can 
build a (1,1) tensor 

XS = sTg™ (2.1) 

that constitutes the basis for all the necessary non trivial scalars [19] 

r n = Tr(X n ) = (XX , n= 1,2,3,4. (2.2) 

Around the background g^ = g^, the combinations (2.2) are clearly invariant under the (global) 
Lorentz transformations defined by g^ v . Because of this residual symmetry, the theories built 
from these invariants were named "Lorentz-invariant" . 3 

For instance, with the popular choice g^ = rj^ and by expanding g„ u = i)u, v + hau on e has 

a (n - 4) 2 + b (r 2 - 2n + 4) = (a h^h^ + bh 2 )+ 0(h 3 ) , (2.3) 
i.e. the Lorentz-invariant generalized Fierz-Pauli mass term. 

We wish to stress that, due to the presence of the fixed external metric g^v, Lorentz-invariance of 
these theories is only present around the particular background g^ = g^. In any other background 
(for instance cosmological) this residual Lorentz invariance is not present. A breaking will appear 
through the interactions in (2.2). In fact, the residual Lorentz symmetry should be not confused 
with the same symmetry entering in the formulation of the Einstein equivalence principle. Indeed, in 
general backgrounds it is not possible to find a locally free fall reference frame such that both g^ and 
g^v coincide with the Minkowski metric. The best one can do is to choose a suitable tetrad basis e° 
such that 

g ab = Diag(-1, 1, 1, 1) = r) ab , g ab = Diag(A , Ai, A 2 , A 3 ) . (2.4) 

Thus, in general, local Lorentz invariance is never present in generic backgrounds, even in the so called 
"Lorentz-invariant" theory. 

Still, if the matter sector is minimally coupled only to g^, in this sector local Lorentz invariance 
is always present, regardless of the invariance properties of the gravitational sector. No Lorentz- 
breaking or violations of the weak-equivalence principle will be observable from the matter sector, and 
any breaking will be present only in the gravitational one, i.e. will manifest itself in propagation of 
gravitational waves, etc. 

In view of this, considering that there are so far no constraints on the Lorentz-invariance of the 
pure gravitational sector, it makes sense to consider also the possibility of a smaller residual symmetry. 



2 General covariance can be also left unbroken by the use of appropriate Stiickelberg fields. In section 6 we will discuss 
in detail this approach, which does not change the result of this section and of our work. 

3 One could also give a dynamical character to the extra metric adding its Einstein-Hilbert action, see for instance [15, 
19, 21-25], leading to a bimetric theory. In the bimetric formulation there is no absolute object but clearly the dynamics 
is more complicated; for instance, at the linearized level there are two gravitons, one massive and one massless, and the 
interaction terms break the "relative" diffs acting on the two metrics. 
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• An immediate simple option is that of a Galileian S0(3) structure, which is achieved by intro- 
ducing in the Lagrangian rotational invariants only. In practice one can build SO (3) invariants 
using the metric components goo, goi and gij (choosing for instance a canonical time-slicing 
x° = const, but see later for comments on dependence on the slicing) or also better using 
the Arnowitt-Dcser-Misner (ADM) [26] variables: lapse N, shift N l , and spatial 3-metric 7^ 



Again, to give mass to the graviton one needs to break spatial diffcomorphisms, so that an 
external spatial metric has to be introduced. For instance, one can use the euclidean cartesian 
metric 5ij. 

As before, unbroken rotations will be a feature of the theory only around the background g^ oc 
Sij. This is however sufficiently wide to include the homogeneous CMB reference background, 
any homogeneous representation of the FRW cosmologies, etc. 

Clearly, this approach includes the above Lorentz-invariant as a subcase, for one can decompose 
the r n under rotations, e.g. T\ = —g + g lJ &ij (for 

9y,v — while in general one has also 
external fixed "shifts" in the game) . It should also be clear that a generic rotationally invariant 
function of the ADM variables can be Lorentz invariant if and only if the variables appear in the 
particular combinations t„. 

Summing up, a general action of massive modified gravity theory can be written as 



Matter is described by L ma ttor and is minimally coupled to the spacetime metric g^, in agreement 
with the equivalence principle. The potential V, taken as a scalar function of the ADM variables plus 
a fixed spatial metric, parametrizes the most general rotationally invariant theory of massive gravity, 
and boils down to the Lorentz-invariant case if V is actually a function of the r„ . 

The analysis in this work will be concerned with determining the allowed form of V to have a given 
number of DoF at nonperturbative level. As noted in [17], the possible Lorentz breaking nature of 
the resulting potentials opens a new window in the space of modified gravity theories to overcome the 
critical drawbacks of the Lorentz-invariant case (most notably the dramatic strongly-coupled nature of 
the theory at phcnomcnologically interesting scales). The necessity of breaking Lorentz symmetry in 
the gravitational sector to have a healthy large-distance modification of gravity should be considered 
as an interesting prediction to be experimentally tested. 

It is worth stressing that with the non-dynamical metrics, which are effectively aether-like absolute 
objects, preferred frame effects must be expected. Indeed, the frame where the frozen metric assumes 
the form 77^ is special and must be specified in relation with matter frame. For instance, one could 
take such a frame as the one where the universe is homogeneous at large scale (the CMB frame); this 
choice is particularly natural if massive gravity has something to do with dark energy. Once such 
a physical choice is made, one should compare against standard preferred frame constraints at the 
post-Newtonian level [20]. In the existing literature those effects apparently have been over looked. 



(U=l,2,3) 




(2.5) 




(2.6) 
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3 Canonical Analysis 



Let us review the canonical analysis for a generic massive gravity theory given in [3] . We use the ADM 
Hamiltonian formulation in terms of lapse N, the shift N l and the spatial 3-mctric 7^ . The presence 
of the deforming potential violates diff invariance and thus the canonical analysis drastically differs 
from the one in GR. Following the very same steps of the ADM analysis, the Hamiltonian relative to 
the gravitational part of action (2.6) is 

H = M| L J d 3 x [N A Ha + m 2 Ny^f V] , (3.1) 

where we introduced the notation N A — Na — (N, N*) and U A = (U, Hi), with A = 0,1, 2, 3. The 
phase space is a 20 dimensional space of the conjugate variables (7^, IF-?) and (N a ,Ua)- We will 
denote by F : a the derivative dF j dN A of any function F, and we define V = m 2 N ^/j V. As matter of 
fact, Ha are functions of the spatial metric 7^ and its conjugate momenta fi' 3 only. Their (standard) 
form is given for completeness in appendix B. 

When m = 0, the action (2.6) reduces to the familiar ADM GR Hamiltonian derived from the 
Einstcin-Hilbert action. In general no residual symmetry of potential V will be assumed, thus our 
analysis applies both to the Lorentz invariant (LI) case as to any nonlinear version of Lorentz breaking 
(LB) models considered in [12, 13]. 

Peculiar of GR and of also the deformed version (3.1), is that TV are non dynamical and their con- 
jugate momenta Ha are vanishing on the physical phase space (the are called primary constraints). 
The total Hamiltonian is 

H T = H+ [ d^xX^A, (3.2) 



where a set of Lagrange multipliers X A arc introduced to enforce the primary constraints « 0. 4 
The time evolution of any function defined on the phase space is given by: 

={f(t,x),H T (t)} = {f(t,x),H(t)} + I d 3 y\ A (t,y){F(t,x),U A (t,y)}. (3.3) 



dt 

The Poisson bracket between two functionals F and G of canonical variables q a , p a is defined as 

5F SG 6F 6G 



{F, G}=J d 3 y 



Sq a {t, y) Sp a (t, y) Sp a {t, y) 8q a {t, y) 



In the following the dependence on space and time of the various fields will be understood. 

Before studying the equations of motion of the dynamical variables hij, it is crucial to enforce 
that the four primary constraints are conserved in time. This leads to secondary constraints Sa'- 

Sa = {Ha,H t } = -(Ha + V,a)^0. (3.4) 

Basically, these are the equations of motion of TV . They do not determine any of the Lagrange 
multipliers and are not automatically satisfied. Therefore, Sa ~ reduce the dimension of the physical 
phase space by 4. Also <!?4 must be conserved, leading to new conditions, the "tertiary" constraints 

Ta = {S a ,H t } = {S a ,H}-VabX b ^0, with Vab = V, A b = q^ n b ■ (3-5) 



The customary notation II j\ f» means that II4 are only weakly zero, e.g. zero on the constraints surface only. 
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If the Hessian |Vab| is non-degenerate, it can be inverted and the four eqs. (3.5) can be solved for all 
the X A ; all constraints are consistent with time evolution and the procedure stops. 

Summarizing, in the general case for which det \Vab\ 7^ 0, we have 4 (TL A )+4 (Sa) — 8 constraints, 
for a total of (20 — 8)/2 = 6 propagating DoF. Technically, these 8 constraints are second class, being 
Rank I {1X4, Sb}\ — Rank|V,4_B| = 4. As a result, the whole gauge invariances of GR are broken. 

In the LI case Lorentz invariance tells us that the propagating DoF must be grouped in a massive 
spin two (5 DoF) representation plus a scalar (1 DoF). This is the Boulwarc-Deser result, valid for 
a generic potential. The extra scalar, the so called Boulware-Deser sixth mode [2], turns out to be a 
ghost rendering these generic nonlinear theories unviable. 

Thus, the first important result is that a necessary condition to have less than six propagating 
DoF, is that r = Rank|VAs| < 4. 

3.1 The case r = 3: five DoF 

In this important case the matrix Vab has a null eigenvector x A an( i three other eigenvectors E A 
(n = 1, 2, 3) with eigenvalues K n : 

Vab X B = , V AB E* = n n E A . (3.6) 

For instance if det(Vf,) 7^ 0, then \ A = (1,—V^ Voj)- The Lagrange multipliers can also be split 
along \ A an d its orthogonal complement 

\ A =zX A + i2 d » < = ZX A + \ A - (3.7) 

n—l 

The Hessian is now non-degenerate only in the subspace orthogonal to \ A S an d from (3.5) one can 
determine only three out of four Lagrange multipliers, d„ — E a {Sa, H} / ' n n . The projection along 
\ A leaves z undetermined and (3.5) gives a new non-trivial (tertiary) constraint 

T= X A {Sa,H}^Q. (3.8) 

Here and in the following we suppose that the constraints found are non trivial. At this point, either 
the conservation of T allows to determine the remaining multiplier z or it may generate a further 
constraint. One has the "quaternary" constraint 

Q = {T,H T } = {T,H} + {T,\ A -U A }, with X A U A = J d 3 y X A (y)U A (y) . (3.9) 

Clearly, the first piece of (3.9) does not contain z. Using (3.6) and (3.7) after some computation (see 
C for more details) we find for the second piece 

{T, X A ■ U A } » x B {\ A ■ S A , S b } - X B {X A Vba, H} + X D V BD - & • z , 

Q(x, y) = X A (x) {S A (x),S B (y)} x B (y) ■ (3.10) 

It turns out that, see appendix B, @(x, y) = A l (x,y)diS^ (x — y) with A(x, y) = A(y,x); it is worth 
to point out that in a finite dimensional context A 1 would be absent. Thus, 



Qz= d 3 y e(x,y)z(y) = --—d x [z(x) 2 A l (x,x)]. (3.11) 

2z{x) 
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As a result, Q is free from z if A l (x, x) = 0, which consists in the following condition 

/v, + 2 x Vp=0, V = 7 1/2 V. (3.12) 

If this condition is satisfied, Q in (3.9) is a quaternary constraint, which together with T confirms the 
elimination of one propagating degrees of freedom (DoF). If the condition (3.12) is not satisfied, (3.9) 
is not a constraint but can be used to determine the last Lagrange multiplier z. This leaves us with 
an odd number of second class constraints and an odd dimensional physical phase space, in particular 
5 + 1/2 DoF. A similar situation is encountered in other deformations of GR, see for instance [27]. It 
not clear if such a peculiar occurrence is physically acceptable or not. For sure, condition (3.12) is 
rather restrictive. 

When (3.12) is satisfied by V, the time evolution of the new constraint Q is non trivial and can 
finally determine the multiplier 

z = -{ X A Q,a)~ 1 [X B Q,b+{Q,H}] . (3.13) 

In conclusion, for r — 3, there arc 4 (ILa)+4 (Sa)+1 (T) + l (Q) = 10 constraints, giving (20— 10)/2 = 5 
DoF, that is a good candidate for a theory of a massive spin-2. The required conditions are (3.12) and 

Rank|V AB | = 3 and x A Q,A ^ ■ (3.14) 

To our knowledge, the observation that a singular Hessian is a necessary condition to avoid the 
Boulware-Deser argument which leads to six propagating DoF was made for the first time in general 
terms in [7]. We would like to note that, though it is not straightforward, the same conclusion is 
reached if we allow for an explicit time dependence in the potential V. 

3.2 General case 

Let us briefly discuss the generic case in which the Hessian has any rank r between zero and four. In 
this case Vab has 4 — r null eigenvectors Xa > 

V AB Xa=0, a = l,...,4-r, (3.15) 
and as before, we can split the Lagrange multipliers as 

4— r r 

\ A = J2 Z <*X A + Y, d " E n- ( 3 - 16 ) 
a—1 n—1 

The projections of (3.5) along E^ determine the r Lagrange multipliers d n = E^{S A , H} / K n , while 
the projections along the Xa gi ye 4 — r tertiary constraints 

T a = X A {S A ,H}^Q. (3.17) 

The conservation of these constraints leads to 

r 4 — r 

Qa = {r a ,H} + J2{Ta,n A }-d n E^-J2^P ■ Zp (3.18) 

n=l $=1 

which are 4 — r relations linear in the remaining 4 — r Lagrange multipliers z a . The number of DoF 
thus depends on how many multipliers are determined by these conditions. The matrix 9 a p is 

e a0 (x,y)= x A (t,x){S A (t,x),S B (t,y)} X p(t,y). (3.19) 
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Notice that due to the dependence on x, y, the matrix 9 a p is not necessarily antisymmetric and its 
rank s is not always even. If s < 4 — r, some z a are undetermined and we have 4 — r — s new quartic 
constraints which reduce the DoF. 

At the next (fifth) step, if the conservation of quartic constraints determines all the needed mul- 
tipliers, the procedure stops and we have only 16 — 2r — s constraints: 4 (Ha) + 4 (Sa) + (4 — r) (T a ) + 
(4 — r — s) (Qa). This implies that the maximal number of DoF is 

DoF<4 + r/2. (3.20) 

This is the maximal number because if some of the z a are not determined, further steps are necessary 
possibly reducing the number of DoF. In general also first class constraint may be present corresponding 
to residual gauge invariances, reducing even more the number of DoF. 

Summarizing, massive deformations with 5 DoF exist only in two cases: r = 3, s = and r = 2, 
s = 2. In the case r = 3, the condition s = is equivalent to the differential equation (3.12) 
for the deforming potential. In the case r = 2 and s = 2, (3.18) determines the remaining two 
Lagrange multipliers and all 10 constraints are consistent with the time evolution. Notice that it is 
not possible to build a potential with s = 2 without breaking spatial rotations. Thus, we have shown 
that r = Rank(V/is) = 3, together with (3.12) are necessary and sufficient conditions for a rotational 
invariant massive gravity theory with 5 DoF (again we consider that none of the constraints found is 
accidentally trivial, i.e. ps 0). 

4 Solution for five DoF: the Monge-Ampere problem 

We have shown that the potentials that have 5 DoF at nonperturbative level are solutions of the 
following system of coupled differential equations in the space of fields (N A 7 7^) 

det(V AB ) = 0, (4.1) 

Vi + 2£ A e'^=0, (4.2) 

where we defined £ A = x A /x° and V = 7 _1 / 2 V. Since spatial rotations are preserved, det(Vy) 7^ 
and the r = 3 condition can be rewritten into the equivalent form 

Voo - VoiCVy)-^ = . (4.3) 

Physically, this means that once the equation for the shifts N l are solved and inserted back in the 
potential, V will be linear in N. Mathematically, (4.1) is an homogeneous Mongc- Ampere equation [28] 
on the space of the N A variables, while jij is a spectator variable. Particular solutions of this equation 
arc known in a closed form and, remarkably, also the general solutions can be given in parametric form, 
which we wil exploit below. Thus, there exist a large class of potentials, taken as a function of lapse, 
shift and spatial metric V (N , N l , 7^) that satisfies (4.1). Note also that if V satisfies the Monge- 
Ampere equation also V does, so we can substitute V — > V into (4.1). Let us describe the general 
parametric solution. 

Our starting point is the existence of a null eigenvector \ A such that \ A Vab = 0. Then it is 
crucial to change variables from N A to \ A L > an d study the properties of the function 

U( X A ) ^VaX A , (4.4) 
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We easily see that U is an homogeneous function of weight one, i.e. 

^K=Va -+ U = x A ^ -+ U(«/)=«U(/). (4.5) 
Therefore we can write 

xj(x A )^ x °m i ), c = K- ( 4 - 6 ) 

X 

where we have introduced a generic function U(^ 1 ). Everything can now be expressed in terms of U 
and relations (4.4), (4.5) become 

which allow one to reconstruct V for any given U. Note, in what follows we will denote with Ui, Uij, 
etc., the corresponding derivatives with respect to the £ l (while instead Vi = dV/dN l ). 

These equations have to be integrable, i.e. the following conditions have to be imposed: 

Vij = Vn and V Qt = V lQ . (4.8) 

The first of these equations boils down to 

which arc known equations describing a generalized hydrodynamic flow. They can be solved formally 
by means of an implicit relation between the N l and £ J , which is, using also the second equation (4.8), 

N l = NC + Q\ Q\e ,7 ij ) = - (U l3 y l £j (4.10) 

where £(£ l , 7* J ') is a new arbitrary function (with again £i = d^£), while (Uij) 1 is the matrix inverse 
of the hessian of U. This is a crucially important change of variables, and should be thought as defining 
implicitly £ l (N l ). 

Intcgrability being assured, one can find a potential V by integrating (4.7). Remarkably, the result 
of the integration can be performed explicitly (see Appendix D) and the result for the potential is 5 

V = NU+Ui Q i +£. (4.11) 

or 

V = U + N- 1 (Ui Q l +£) . (4.12) 

This completes the general solution of the Mongc- Ampere equation. 

Summarizing, the solutions can be parametrized with two arbitrary functions U and £. Concretely, 
the strategy to generate explicit potentials function of (N, N l , -f lJ ) is therefore first to choose U and 
£ which already fully determine the potential, and then to solve (4.10) for £ l in terms of (N l , N, Jij). 

As a final step, we have to impose condition (4.2) in order to fully reduce to five DoF in phase 
space. This brings j lJ into the game. It is nice that, when expressed in terms of ZY, (4.2) takes an 
elegant and compact form 



5 Note: during the integration of (4.7), the choice of an integration constant C(7 U ') is required. This can be reabsorbed 
in £ without changing the previous steps. With respect to [17], we have here £ = C — C. 
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whose solution is readily found as 

U = U(lC ij ), with JC ij = j ij - . (4.14) 

This completes the solution of our problem, which is described in terms of the two functions U (/C y ) 
and £ 

Since we have assumed unbroken rotations, 1A must depend only on scalar combinations built with 
JQ 3 and (and not with 7^ because this would invalidate cq. (4.13)). So, the rotationally invariant 
solutions of eq. (4.13) are of the form U (Tr[(/C/)] , Tr[(/C/) 2 ], Tr[(/C/) 3 ]). 

For generic U and £, the resulting potential (4.12) will in general be Lorentz breaking. It is in fact 
clear the reason why a full Lorentz symmetry is not natural: the basic scalars that solve (4.13) are just 
three, while there are four SO(3, 1) invariants built out of X£ = g^ ia rj au , for instance its eigenvalues. 
Thus, in building the potential V from hi, the four 50(3, 1) scalars should come out only with a very 
special choice of the function £ entering (4.10) (e.g. see (5.15) below). This is in agreement with the 
discussion above, in section 2. 

Before describing explicit examples we can also evaluate the hamiltonian H, after the use of the 
(secondary) constraints (we neglect the surface term that is formally identical to the one of GR): 



H constramed = J d 3 x (n A N A + V) SA ^ ^ j d 3 x (-V A 2V A + V) eq(47) 

d 3 x 1 1 ' 2 (y-uN-UiOfy 

= Jd 3 xj 1/2 £. (4.15) 

Thanks to this remarkably simple expression, one can have positive definite energy by choosing ap- 
propriately £ . 

5 Explicit theories with five DoF 

Now that we have parametric solutions to the constrained Monge- Ampere equations as free functions 
of U and £, we have to Seattle last step with the explicit calculation of the potential V through 
eqs.(4.7). Let us look for some specific examples where the potentials can be given in a close form. 

5.1 Lorentz breaking solvable examples 

Case 1. A simple but quite wide set of theories is found by looking for £ such that Q 1 cx £ 4 : 

Q l = Ch)€ (5.1) 

for any arbitrary scalar function £ of > so that 

N* .... ATW- 7 ' 
P = and V' J = Y J - 7 ^ ■ (5-2) 

This is realized by the following £ 

£ = - (£% -U) C + C, (5.3) 
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for any function C of 7 y . From this we find the potential 

V ={N + QU{K)+C{ 1 ). (5.4) 

The minimal case of ( = was considered in [17] and it already has a very interesting phenomenol- 
ogy. We noted there that the choice ( = together with C = is undesirable as it leads to strong 
coupling in the scalar and vector sector. 

Case 2. An other solvable case is if Q l is linear in 

Q l = M) (3 (5.5) 

with M a matrix that can depend on 7**7^-. Thus, (4.10) can be inverted to express £ in terms of the 
lapse and shifts 

ff= (N + M)- 1 ^. (5.6) 

where it is here convenient to use the matricial forms N = {N 1 }, £ = as well as U" = {Uij}, 

£' = {Si}, etc. For £, eq. (5.5) implies that 

£' = -U" M£, (5.7) 

whose integrability conditions £y = arc 

U" M — M* U"=0. (5.8) 

If M is not the identity (which case corresponds exactly to case 1 above) this condition can be satisfied 
only if U" is independent of £. Thus U must be at most quadratic in £ l , and because it must be a 
function of JC % ' J = 7 IJ — ^'f , the only allowed term is Tr(£ • /). As a result we have 

U = C - £ Tr(/C/) = C - ^ TV(7-V) + % ?f?; 

J 2 2 (5Q) 

where C, C are functions of 7 and c\ is a constant. The potential V can be derived for general M along 
the lines of the previous examples and clearly gives an expression which is quadratic in the shifts 



N T (N + M)- 1 (f + N- 1 M-f)(N + M)- 1 N-Tr(-/- 1 f) +C + N- 1 C. (5.10) 
The corresponding on-shell Hamiltonian (4.15) reads 



H = j d 3 xj 1/2 — y N ■ (N + M)~ l ■ f ■ M ■ (N + M) _1 ■ N + C . (5.11) 
5.2 Lorentz Invariant case 

Let us consider now the Lorentz invariant case. The change of variables required for the solution of 
the Monge Ampere, 

N' 1 = NC + Q l (tl) (5.12) 

coincides in fact with the change of variables used in [6] to show that the Lorentz invariant potential 
found in [5] propagates five degrees of freedom, (4. 10), after relating the variables n l to £*, via « l (£, 7) = 
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(D- 1 )*^'. This gives a nontrivial implicit relation between N l and £*. In detail we have (we limit 
ourselves to / M „ = 77^) 



2\ ' 



[l + e-zc- 1 -^ (/c-/)j ; e i = (x>- 1 -o < 



with £ = f J - f f and AT 1 = 7^ 



i-e 



and where £j = 7^ £ J and £ 2 = £ l 7^- £•? . The two unknown functions U and £ are 



W L7 = Tr 



(£•/) 



Computing = — [(/C • /) ' ■ we have that 



1/2 



3, £ = (i+e-/c- l .c) 1/2 = 



and the full potential turns out to be 



1/2 



V = N Tr 



(£■/) 



1/2 



3) +^1-^ = V 7 !"^ (l + iVTr(£>)) -3JV. 



(5.13) 
(5.14) 

(5.15) 

(5.16) 
(5.17) 



Here we see that the condition £ 2 < 1 must hold to have a real action. 

Going back [6] to the original variables, we find the Lorcntz invariant expression in terms of 
XH = Q^fau, namely 

Vli = MX 1 / 2 ) - 3 , (5.18) 

Thus, clearly also in the Lorentz invariant case, U is a function of JC only and it satisfies automatically 
the nontrivial condition (4.13) for the complete elimination of the sixth mode in phase space. The 
on-shcll Hamiltonian is 



H 



(5.19) 



and it is real and positive definite only for £ 2 < 1. Notice that actually when £ 2 = 1, K % i is not 
invertible and the change of variables from £ l to N 1 is ill defined. 



6 Stiickelberg Dictionary 

Even if the theories considered in this work feature a massive graviton and thus a direct violation 
of general covariance, this invariance can be restored, or considered never broken, by the use of 
Stiickelberg fields. This can be achieved by promoting the external metric(s) or tensorial quantities 
to true tensors, in such a way that the action is fully cliff invariant. 

For instance, as is well known, in the case of Lorentz invariance, one can promote the external 
metric in (2.1) to a true tensor by using four Stiickelberg fields $> A (A = 1, . . . , 4): 

9nu = d^ A d v <\> B tub , (6.1) 

after which all traces of X£ are also diff invariants. 

To promote rotational scalars to full diff invariant quantities one can introduce a similar set of 
Stiickelberg scalar fields, starting with "universal time" field $ which defines a spatial slice through 
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<3> = const, with unit normal = g a Pd a $dp$) 1 ' 2 . In the spatial hypersurface $ = const. 

one also considers the 3D metric that is pulled back to spacctime with three more Stiickclbcrg 
scalars $' 

I,, f ij d lt &d v &. (6-2) 

Of course, the presence of a residual SO(3) symmetry in the unitary gauge requires that = 5ij. 
Note that in this procedure, the preferred time slicing has become dynamic with the gauge field <E>. In 
the unitary gauge, x° = $ and x % = $ z we have that 

g^ = 6 ^ d^ = S l = (6-3) 

In general, any 50(3) scalar in the unitary gauge can be written out of the following 4D basic 
scalar, vector and tensor quantities 

n=(-^d a Zd^)-V\ ^—A*—^, Yf-gTU. (6.4) 

I.e. the r n of (2.2) can be decomposed in terms of new scalar combinations of them, for instance we 
have n = ?i~ 2 + g^ffiu- 

After this dictionary, all the potentials discussed in the previous sections can be easily written 
in a generic coordinate system simply by noting the the lapse N should be replaced by n, all scalar 
functions of Y k $kj replaced by the same scalar functions of T£ defined as 

r£ = (g^ + n a ) f av = Y<? + n" n a f a „ . (6.5) 

which plays the role of the pull-back of the 3D spatial metric. 
The potential for the case 1 is given by 

l/ caao l = n- 1 [n + Q{T)]U(fc) + n _1 C(r) . (6.6) 

where 14 is now a scalar function of 

^ = ( r "° ~ V+Wf) / °" ' (6 ' 7) 

The expression for V in the case 2 is not particular illuminating and it can be determined along the 
above lines. 

Thanks to the presence of the Stiickelberg fields, the actions maintain the invariance under dif- 
feomorphism. Once the unitary gauge is chosen, only a residual Lorentz or rotational invariance is 
left. 



7 Perturbations around flat space 

From the weak field expansion of the action (2.6), taking g^ v = n^ + h^, at the quadratic order in 
h we have 



S = I d 4 xM^ 



£(2) + — ( m o h oo + 2to i h oihoi - ml h l3 hij + mj h\ - 2m\ h 00 h u ) 



(7.1) 



where £n) is the standard Lagrangian for a massless spin 2 particle in Minkowski space, stemming 
from the Einstcin-Hilbert action expansion. The other terms represent the most general form of the 
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expansion of any potential with a SO(3) residual symmetry. The physical consequences of a such mass 
term were first discussed in [12]. 

Here, using (4.12) the masses can be computed explicitly by differentiating the potential. As 
the first step, one has to find the conditions for the potential to admit the Minkowski background 
SV = Viif These are 



U = 0, 



(7.2) 



The bar indicates that all expressions are evaluated on Minkowski space and dS/d^ = £'jij, 
dU/dJC J = U'jij. We have also used the fact that by rotational invariance on Minkowski space 
Qi\ v = Ui\. n = C\ n = 0. From the general expression of the potential (4.12) and the condition (7.2), 



we find the masses 



m 2 d 2 (^gV) 



4 ON 2 
m 2 d 2 (^gV) 



= —6i 



dWdN 1 



m 2 U' 



( d? 



(7.3) 



djvdN 



3m< 



w 



The expressions for m| and m§ are not particularly illuminating. It is reassuring that the condition 
m\ = 0, to have only 5 DoF [2] at the linearized level, is automatically satisfied. Moreover, the 



following relation holds 



dN* 



(7.4) 



In a Lorcntz invariant theory where ra\ = vn\ one should have that (d^/dN l ^\ = 3/2, as can be 
checked from section 5.2. In the Lorentz breaking case, a large class of potentials easily passes the 
solar system test and are weakly coupled in all phenomenological interesting scales [17]. 

Let us briefly review the linearized analysis of [12] (see also [13]). It is convenient to expand the 
perturbation in terms of 50(3) scalars, vectors and tensors, namely 



h 



00 



div, 



Xij 



didj 



djQi 



didj a 



(7.5) 



It is easy to see that the transverse and traceless tensor \ nas always a standard kinetic term and a 
mass m|; thus the only requirement is m 2 , > 0. In the scalar sector, when m 2 , = 0, ip is non-dynamical 
at appears linearly in the quadratic action. Again, from (7.3), when mi ^ also ^ and we can 
solve for a. We are left with a single propagating scalar. In the vector sector m is non dynamically 
and can be integrated out, leaving a single propagating vector a^. As a result, we have five DoFs that 
match canonical analysis. On the other hand, if mi = 0, in the vector sector there is no propagating 
DoF. From nonperturbative canonical analysis we already know that actually there are 5 DoFs; thus 
the vanishing of the kinetic term in the vector sector when mi = it is just an accident of the linearized 
expansion around Minkowski space that signals that there is a vector mode that is strongly coupled. 

For a static source with a non-zero energy momentum tensor component too, the gauge invariant 
potential $ = hoo is [30] 



-3 m„ 



$ = 2t 00 - 



4 A 2 



6 m| 



^oo 
2 



A - M\ 



M 2 



(7.6) 
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where the squared masses Mf 2 and the coefficients A\ 2 in the formal decomposition of the last term 
depend on the mass parameters mf , with A\ + A2 = 1. It is straightforward to check that under the 
mild conditions 

0<ml<lml, "2^2 + ! y + ! y <m 2 ( 7 - 7 ) 

the squared masses Mf and Mf are real positive, so that the potential is a sum of two Yukawa terms: 

§ = -(A! e" Ml r + A 2 e~ M * r ) (7.8) 
r 

Also, if one of the A\^ vanishes, the potential mimics a single Yukawa term of the FP theory. 

The conditions (7.7), via (7.3) and the similar expressions for m2, m,3, turn into point-wise con- 
straints on the derivative of V on Minkowski space, which can be always satisfied by appropriate 
choices of the two otherwise arbitrary functions IA and £ . 

By (7.6) and by inspection of the other gauge invariant potential in the metric fluctuation, one finds 
that both potentials reduce to 1/2 A at short distance, confirming the absence of vDVz discontinuity 
[12]. The phenomenological advantages of having a weak coupled theory compared with the Lorentz 
invariant case has been discussed in [17]. 

8 Conclusions 

In this this paper, by using the canonical Hamiltonian analysis, we studied the number of propagating 
degrees of freedom in generic theories of massive gravity. Such theories which lead to large distance 
modifications of gravity, include nonderivative terms in the action, encoded in a potential V function 
of the metric components. To build V extra tensorial objects are needed. We show that, in order 
to preserve at least the rotational symmetry SO(3) and lead to a massive graviton, at least an extra 
3-dimensional metric is needed. In the minimal approach such a metric is non dynamical, and will 
lead to potentially observable preferred frame effects. 

From the canonical analysis, it follows that 6 DoF propagate in general unless additional conditions 
on V are introduced. In particular we have shown that 5 DoF propagate if and only if the potential 
satisfies the Mongc-Ampcre equation (4.1) plus the additional differential equation (4.2). Such set of 
differential equations can be solved in terms of two arbitrary functions which parametrize the most 
general nonderivative gravity modification with 5 DoF at full non perturbativc level. Thus, by using 
the canonical analysis, the problem of finding the most general theory of modified massive gravity 
with 5 DoF has been completely solved. Let us summarize the result in a nutshell: given two arbitrary 
functions W(AC* J ') and £(£*,7* J ') (of the specific arguments) the deformation potential V of any massive 
gravity theory with five propagating DoFs can be written as 

V = U + N- 1 (£+Ui Q l ) where 

N* = Nc+Q i , q'ix'---' j ) K ij = -f ij - e e . 

In addition, £ represents the energy (Hamiltonian) density of the system and thus it can be chosen to 
be positive definite. 

Besides its theoretical interest, this result is also relevant from a phenomenological point of view, 
for a number of reasons. First, because it uncovers a large class of massive gravity theories that are 
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ghost free on Minkowski space, whereas previously, the only known ghost free theory was the Lorentz- 
invariant four parameter theory of [5, 6], which is shown here to be a very special case of our general 
construction. 

Second, while Lorentz symmetry may be enforced, the price to be paid is the impossibility of using 
perturbation theory in many physical important situations like inside our solar system. Instead, in 
generic Lorentz-breaking theories this problem disappears. It ought to be remarked that the Lorentz 
symmetry we are discussing here only concerns the gravitational sector, and is not the same symmetry 
that enters in the formulation of the Einstein's equivalence principle. As such, is not subject to 
phcnomcnological constraints. Thus, we conclude that the viability of the theory is directly connected 
with the need to have Lorentz-breaking in the gravitational sector, which should also be testable at 
forthcoming gravitational wave experiments. 

The concrete phenomenology of the new class of Lorentz breaking theories is also rather promising, 
as argued in [17]. From a perturbative point of view, from the general expression of V, we have 
found that remarkably some relations among the graviton masses exists and are important in view 
of the phcnomcnological application. From a nonpcrturbative side, besides the absence of ghosts 
in the spectrum, crucially important features are the possibility to trust the theory up to the cutoff 
A 2 = (mM PL yi 2 ~ (lO" 3 mm) , and the absence of strong nonlinearities (Vainshtcin radius) around 
macroscopic sources. 

It is crucial to study the impact of the general construction presented here on cosmology. We 
leave the matter for a future publication [31]. 

A Canonical analysis made simple in the rank 3 case 

In the Monge- Ampere case, the canonical analysis is made very simple and transparent if the £ l 
variables are used in place of the shifts N l , through the change of variables 

N i = N? + Q i (C,i ii ). (A.l) 

The Hamiltonian, as a function of the (N, £) variables, is 

H = N{H + CUi) + Q l U t + V(N, NC + Q l ) . (A.2) 

The primary constraints are as usual ILa = dH/d(N,^ 1 ) as 0, and one introduces four Lagrange 
multipliers X A and a total hamiltonian Ht = H + A^ILa. 

The secondary constraints are now modified, but after using the Monge- Ampere relations (4.8) 
and (4.10), one has the following simple expressions 

50 = Ho+eHi + y/iUaO; (A.3) 

51 = Wi + V7«i«0. ( A -4) 

In the simplified form the secondary constraints are expressed in terms of the function U, which 
depends only on £ and not on N. Notice that So = Sq + £ l Si as Sq . The tertiary constraints are 

To = {S ,ff}«0; (A.5) 
% = {Si,H} + VijX j w . (A.6) 

The Lagrange multiplier A° drops out from the expressions, and this is due to the fact that U on So 
does not depend on N . Thus, To represents a genuine new constraint, eliminating half of the sixth 
mode in phase space. 
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Finally, the time evolution of 7o is 

Q = {T 0} H T } = {5 , 5 } A + (A°-independent terms) ; (A.7) 

and if one requires that {So, So} = 0' wc S c ^ a further constraint which completes the elimination of 
the sixth mode. This condition, using (A. 3) and the algebra of GR (see appendix B), leads straight- 
forwardly to the additional condition (4.13). 

Thus, the implicit change of variables (4.10) from N l to £ l has the effect of reducing the hessian 
null eigenvector to (1, 0), and it helps to make more transparent the generation of the new constraints. 



B Poisson Brackets 

We calculate here explicitly the poisson brackets required for 

Q(x,y) =x A (x) X B (y) {S A (x) 7 S B (y)} 
=X°(x)x (y){So(x),S (y)} 

+ X°{x) X 3 {y) {S (x), S,(y)} + X l (x) X °(y) {$<(*), S (y)} 

+ x l (x)x J (y){s t (x),s J (y)}. 

Using the following standard expression H and "Hi for the Hamiltonian and momentum constraints in 
GR 

u(x) = -L [2n m „n™ - ( 7 iJ 'n y ) 2 ] - ViR {3) ; 

V7 (B.2) 
Hi(x) = -2 lim D k U mk ■ 

we have that 

{S (x), S Q {y)} = W{x)d\ x) 8^{x-y)-n i {y)dl y) 5^{x-y) 

+ [T mn (x)V^(y) - F mn {y)VT{x)} S i3 \x - y) ; 

{S (x), SM}=ny)d ( f ) ^\x-y) + 6^{x-y)T mn {x)VY' n {y) 

+ V a N b (x) Q ab] (y)5 i3) (x - V) -2V a N \x) -f j{a (y)d^S^(x - y) . 
{Si(x), S (y)} = -H(x)dl v) 6^(x^y) - 6® (x - y) F mn {y)V™ (x) 

- V&(jt) Qabi(x)5^(x -y) + 2 V$(y) li{a {x)d$ 5& {x - y) . 

{S % {x), S 3 {y)} = n ] (x)dt ) 6^(x -y)- Hi(y)a^6^(x - y) 

- 8&(x - y) Vf(y) Q abi (x) + 2Vf(y) ll(a {x)d\* ] 5^ (x - y) (B.6) 
+ 5^(x - y)Vf{x) Q ab3 (y) - 2Vf(y) lj(a (y)d { b f 8^ (x - y) . 



where 



(B.3) 
(B.4) 
(B.5) 



F mn {x) = 7~ 1/2 {2Ii mn {x) - K(x) lmn (x)) ; 

Qabrn(x) = d^-f ab (x) - 2d^J ^ b)m {x) \ /g^x 

Vmn Y \)tnn 



dN8 lmn ' 1 dN^ 
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We have also used the following relation for the Lie derivative C of spatial metric with respect of a 
spatial vector e 

(£ e7 ) a6 - D a e b + D a e a = e m d mlab + lmb d a e m + lam d b e m . (B.8) 
Assembling the various pieces we arrive at the result 

B(x, y) = A\x, y) d\ v) 5^\x -y) + B\x, y) d\ x) 8^{x - y) + C(x, y) 6<® (x-y); (B.9) 

with 

A k (x, y) = - x° (x) X ° (y) n k (y) - 2 X ° (x) X 3 (y) l ja (y) V a N k (x) - X l (x) X k (y) U { (y) 

- 2 X \x)x\y)j 3 a(y)v? k (y) - x k (x)x°(y)n(x) ; 

B k (x, y) =x\x)x°{y) H k {x) + X k (y)x°(x) U{y) + 2 X (y) X l (x) lia (x) V a N k (y) 
+ 2 x k (x)x 3 (y)H 3 (x) + x l (x)x 3 {y) lia{x) Vf(x) ■ 



(B.10) 
(B.ll) 



and 



C(x, y ) =x°(x)x°(y) [Cmn(x)VT{y)-Cmn{y)VT{x)\ 

+ X°(x)x l (y) [C mn (x)Vr(y) + V a N \x) Q abl (y)} 
~ X\x) x °{y) [C mn (y)Vr (x) + V$(y) Q M (x)] 
+ X\x)x 3 {y) [Vf{x) Q abj (y) - Vf(y) Q abi (x)] 



(B.12) 



Thus 



I(x) = J d 3 yz(y)e(x,y) = D i (x,x)d i z(x)+z(x)d$ y) D i (x,y) ]y = x ; 
D i (x,y)=A i (x,y)-B i (x,y). 

The same result can be also obtained by using d^S^(x-y) = -d^S^(x-y). 
We notice that 

C(x,y) = -C(y,x) => C(x,x) =0; 

B\x,y) = -A\y,x) D l (x,y) = A l (x lV ) + A* fax) => D\x,y) = D\y,x) . 
From the symmetry of D l , it follows that D l (x,y)\ y=x = 1/2 d[ D l (x, x) and as a result 



(B.13) 



(B.14) 



I(x) = J d 3 y z(y)Q(x,y) = ^-^[zix) 2 D\x,x)] . (B.15) 

The integral I{x) vanishes for all z if D l (x, x) = 2A l (x, x) = 0, thus the above condition requires 

2 X A X m V l Aln m ~ X a \f° - X 1 X A Va = , (B.16) 

which is precisely (3.12). 

C Quaternary constraint 

In the case r = Rank^V^s) = 3, the only tertiary T constraint reads 

T = j t ix A s A ) = X A [{s A , n} - x B v AB ] « o . (c.i) 
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We have to impose that also T is conserved 

dT 
dt 



= {T, H} + J d 3 y {T(x), X A (y)U A (y)} . 



(C.2) 



The first two terms do not depend on z; using Jacobi identity, the last term can be written as 



{T(x), \ A (y)U A (y)} = y 



A {S B ,H} + {V AB ,H} X B 
d 3 y\ A (y)x(x) B {S A (y),S B (x)} 



(C.3) 



Now using that A = z X + A , the properties of X an d the previous constraints we have that 



dx 1 



dx 1 



(v BC \ A \ c ) ; 



Q N A >. ' J djV^ 4 



(C.4) 



Thus 



3 j/ {f(x),X A (y)U A (y)} 



dN A 



V BC \ a A c +\ A {V AB , H} x l 



d i y\ A {y)x B {x){S A {y) 1 S B {x)}~ / d s y z(y) Q(x,y) 



(C.5) 



As a result the quaternary constraint does not depend on z if (3.12) holds, and it gives a genuine new 
constraint. 



D General solution of the Monge- Ampere Equation 

The solution of the Monge- Ampere equation V is given implicitly in terms of the two functions 14, C 
of £ and 7 by 

W = U-™P d^ = ou 

ON dp* ' ON' dp' K ' ' 

The relation between N l and £ reads 

N l = N C + Q l (D.2) 
Q i (e,7 ij ) = -(U ij )- 1 £ j . (D.3) 

Integrating the first of (D.l) with respect to TV we can exploit the relation that follows from (D.2) 



dU , „8U dQ l dU 

-p = N h — r 

dp* dN dN dp 



Thus 



Now, from the definition of Q 1 we have that 



d£ dp 



dpdN 



dp dp dN dp dN 



d£ 
dN 



(D.4) 



(D.5) 



(D.6) 
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As a result 

V = NU+£+C: (D.7) 

where C is function of N l and 7^. Finally we have to enforce the last of (D.2), this can be done by 
setting 

C = Q%- (D.8) 
Thus we get the final expression: V = NU + £ + Q l Ui. 
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